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We report on a search for the free decay component of 7r~ stopped in water and light materials. 
A non-zero value of this would be an indication of anomalous 77 e contamination to the Ve. and z7 M 
production at stopped-pion neutrino facilities. No free decay component of ir~ was observed in 
water, Beryllium, and Aluminum, for which upper limits were established at 8.2 x 10~ 4 , 3.2 x 10" 3 , 
and 7.7 x 10 -3 , respectively. 

PACS numbers: 14.60.Pq, 29.25.-t, 36.10.-k, 25.80.Hp 
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I. INTRODUCTION 

Flavor suppression in the incident neutrino beam is 
one of the key elements for successful neutrino appear- 
ance experiments. In the case of beam-dump experiments 
using muon decays from stopped pions, the contribu- 
tion to neutrino production from negative pions is ex- 
pected to be significantly suppressed, because negative 
pions are believed to undergo nuclear capture promptly 
after stopping in material, while positive pions freely de- 
cay. The ratio of ir~ /ir + production cross sections for 
a proton beam in the 0.5-1 GeV energy region is about 
1/5 for light nuclei, which further enhances this asym- 
metry. The total neutrino content after the pion decay 
7r + — > [iv^ followed by the muon decay /i + — ■> e + ^ e 77 M 
(tt + — > /i + — > e + decay) does not include 77 e 's in the final 
state. The contribution from decay-in-flight is a few per 
cent for a pion. This provides a good opportunity for 
V e appearance experiments.Q,0,|j|] However, a finite ir~ 
lifetime in the stopping material may result in additional 
V e contamination. 

Normally, after stopping and being trapped in an 
atomic orbit, a negative pion orbiting the nucleus 
promptly de-excites by the Auger process, radiative cas- 
cade decays, or Stark mixing to an S state, where the 
significant overlap of the wave function with the nucleus 
enables the ir~ to be captured by the nucleus. If the ir~ 
is not captured by the nucleus it may decay freely in the 
atomic state. The free decay fraction of 7r _ 's was mea- 
sured to be 8.7 x 10~ 5 Q in Hydrogen and is expected 
to decrease with the atomic number. However, there 
is an exception in the 7r~He atom system^, where 
the 7r~ is trapped in a meta-stable state that has little 
wave-function overlap with the nucleus and the Auger 
processes are suppressed 0- It is therefore important to 
confirm the absence of other mechanisms that may cause 
free decays of 7r~'s. 



Limits for the fraction of the free-decay component 
exist for H 2 and He0,0,0 (data in Ne is only applicable 
to a meta-stable state with a lifetime longer than ~ 1 
ns). This paper reports the results of a search for a free 
decay component of 7r~ in water and light materials 
from analyses of electron yields and proton spectra. 



II. PRINCIPLE 

Figure 1 illustrates the time diagrams involved in pi- 
onic atom formation (left) up to the disappearance of 
the pion by (A) prompt nuclear capture, (b and B) free 
decay with a partial lifetime t w , and (C) delayed pion 
capture with a partial lifetime of T cap t- The processes 
(a) and (c) feed into "7r~ Atom" (right). Since the di- 
agrams for delayed pionic atom formation (left) and de- 
layed nuclear capture (right) are identical, and since the 
two processes cannot be discriminated or separated in 
this analysis, the two processes are considered as a single 
process with prompt capture notation, i.e. R and T capt 
are used in this analysis. The common observed lifetime, 
T b s = 1/(7- + - 1 - ), in the diagrams already assumes 
only delayed capture process. 

A direct method to search for free ir~ decay is to detect 
decay products of pions (B). Detection of muons from 
the decay 7r~ — > [i~Ts in water is not practical because 
of the short muon range in water (~ 1.4 mm). How- 
ever, electrons from the — > fj,~ — » decay chain 
(E e =0 — 53 MeV) can be detected. In this method, 
the background from decay-in-flight (DIF) of a pion to 
a muon, which stops in materials around the target and 
later decays to an electron, must be subtracted using 
Monte Carlo (MC) calculations. The advantages of the 
search for electrons from the decay 7r~ — » /i~ — > er are 
that the measurement is the sum of all free decay com- 
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FIG. 1: Left: Pionic atom formation processes, and Right: 
Atomic de-excitation and capture processes of the pionic 
atom, r and R are the fractions of delayed pionic-atom for- 
mation and delayed capture, respectively. r/ orm and T cap t 
indicate the partial lifetimes for delayed pionic-atom forma- 
tion and delayed capture, respectively. 



portents with different lifetimes and that the sensitivity 
does not strongly depend on the lifetime t ous ; it varies 
only by ~10 % for the lifetime between ns and the pion 
lifetime r v = 26.033 ns0. 

Alternately, free 7r~ decay may be detected by the pres- 
ence of electrons from the decay n~ — > e~v. However, 
because of its small branching ratio, 1.231 x 10~ 4 0, it 
is very difficult to differentiate it from the background 
arising from prompt nuclear capture 7-rays that convert 
in the target to an electron-positron pair. The search for 
7T~ — > e~v is effective for a component with a relatively 
long lifetime. 

The third method (indirect) is to search for delayed 
components (C) in the time spectrum of pion capture 
products, such as protons, 7-rays, and 7r°'s. The ratio of 
delayed ir~ capture f cap t is expressed as 



fcapt — 



(C) 



R- 
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1-R + R- 



(1) 



By eliminating R and T cap t, which are not directly ob- 
servable, the free decay fraction (B) can be deduced, 
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It should be noted that when t oos approaches tv the equa- 
tion approaches to unity, indicating that the capture ex- 
periment is insensitive to the free decay component as 
Tcapt approaches infinity. In order to be sensitive to a 
fast decay component in the presence of the dominant 
prompt capture (the process (A)), the experiment is re- 
quired to have good time resolution. In the search for de- 
layed components in water, this can be achieved by mea- 
suring protons that have high energy losses, therefore, 
with small time walk and spread. Proton measurements 



are not sensitive to ir~ capture by Hydrogen in water, 
since only 7-rays, 7r°'s, and neutrons are produced. The 
majority (> 99 %) of 7r~'s in water |lfl llll | are captured 
by Oxygen atoms. 
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FIG. 2: Experimental setup. TINA and MINA are Nal(Tl) 
detectors described in the text. 



III. EXPERIMENT 

The experiment was performed at the TRIUMF M9A 
channel. Typical ratios of incoming beam particles at 
the experimental target were e:/x:7T=5:l:4 for a positive- 
particle beam and e:/i:7r=88:3:9 for a negative-particle 
beam at P^ = 70 MeV/c with a momentum bite AP/P = 
±4 %. Figure 2 shows the detector set-up. The position 
and the direction of incoming beam particles were mea- 
sured by two sets of X-Y read-out wire chambers (WC1 
and WC2) located 95 cm and 50 cm upstream of the tar- 
get. Two plastic scintillators Bl (76-mm high, 76-mm 
wide, 3.2-mm thick) and B2 (35-mm high, 25-mm wide, 
5-mm thick) were placed 8 cm and 3 cm upstream of 
the target center, respectively. In order to collimate the 
beam (to detect scattered beam by e.g. the WC frames), 
a 305 x 305mm 2 plastic counter with a 76-mm diame- 
ter hole (not shown in the figure) was placed 5 cm up- 
stream of the Bl counter. Purified water for the target 
was held in a 60 x 60 x 8 mm 3 container made of 0.2- 
mm thick aluminum frame and 50-/xm thick beam win- 
dows (both ends). The target was tilted to 45 degrees 
with respect to the beam (some runs were taken at -45 
degrees). A 152xl52-mm 2 veto counter B3 for the de- 
tection of passing-through beam particles was placed 40 
cm downstream of the target. More than 90 % of pions 
stopped in the target at a rate of 1-10 kHz for the it~ 
beam and 30-100 kHz for ir + . 

The decay products were observed by two telescope 
arms at 5 cm from the beam axis at ±90 degrees. Each 
telescope consisted of three 3-6 mm thick plastic scintil- 
lators (T1-T3 or M1-M3), and a 46-cm diameter, 51-cm 
long Nal(Tl) crystal (TINA) or a 36-cm diameter, 36- 
cm long Nal(Tl) (MINA). In each telescope, the farthest 
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scintillator from the target was the largest one that cov- 
ered the front face of the Nal for pile-up detection. The 
solid angle of each telescope was defined by the 102 x 
102-mm 2 Tl or the 76 x 76-mm 2 Ml counter. For tt° de- 
tection, a 5-mm thick 51 x 51-mm 2 lead sheet was placed 
as a converter between the first and the second plastic 
scintillators in the MINA telescope. 
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FIG. 3: Energy in Bl vs beam burst (RF) time with respect 
to B2 time. The time walk effect in the B2 counter has been 
corrected. 
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FIG. 4: Energy loss in T2 (MeV/cm) vs energy in TINA 
(MeV). Electrons below 45 MeV are from n — ► fx — ► e 
decays and around 60 MeV are from electron scattering. 

The trigger required a pion stopping in the target 
(B1-B2-B3) with the presence of an outgoing particle into 
the T2 or M2 counter in the time window from -50 ns 
to 150 ns with respect to the prompt time — this condi- 
tion allowed to accept events containing electrons, 7-rays, 
7r°'s (two 7-rays), protons, and deuterons. Typical trig- 
ger rates were 70 Hz for n~ runs and 150 Hz for n + with 
a deadtimc of 10-30 %. The timing and pulse heights of 
all the scintillation counters, two coordinates from each 
wire chamber, and the relative phase of the 23.06-MHz 
cyclotron RF (incident proton beam bunches) to the B2 
counter time for the time-of-fiight (TOF) measurement of 



a beam particle were recorded. In order to record pile-up 
information, additional wide-gate (250 ns) ADC's were 
employed using signals split from Bl, B3, T3, and M3 
counters for detection of extra charge outside the stan- 
dard gate (50 ns). Pile-up in the Nal counters was mea- 
sured by using two additional ADC's with an early gate 
that closed just before the main pulse, and a late gate 
(by 200 ns) to compare to the normal gate ADC. 

Special runs were taken for calibration and normal- 
ization; runs with a trigger using a pulser driven by a 
pre-scaled cyclotron RF signal, runs requiring only the 
presence of incoming particles (B1-B2), and runs with 
positive beam particles (for normalization of the electron 
acceptance). For tt — > n — > e measurement, a 30-mm 
diameter, 12-mm thick cylindrical water target was used 
without tilting, together with a 25-mm diameter, 5-mm 
thick B2 counter that faced the target with a variable 
gap between z = — 3.5 cm. Also, other candidate 
materials for the beam moderators at a neutrino factory, 
Beryllium (50 x 50 x 8 mm 3 ) and Aluminum (75 x 75 x 6 
mm 3 ), were used as targets during several runs. A 
vinegar target (5 % acetic acid by volume) was used to 
simulate possible effects of excessive H+ ions in water 
generated by the high-intensity proton beam of the 
stopped-beam facility. Measurement was made at beam 
momenta of 64, 67, 70, 73, 76 and 80 MeV/c to confirm 
the stopping position of pions and to test systematic 
effects. 



IV. ANALYSIS 

Beam pions were identified by the energy loss in the 
first beam counter Bl, and the TOF with respect to the 
proton beam bunch. Figure 3 shows a scatter plot of 
Bl counter energy vs RF time with respect to the B2 
time after the time walk correction. At the experimental 
target position (about 10 m from the pion production 
target), the TOF of electrons coming from the next beam 
bucket was separated from that of 70-MeV/c pions, and 
the probability of misidentifying a beam electron as a 
pion was estimated to be - 2 x 10~ 7 . The TOF peak of 
muons was clearly separated. 

The beam counters were calibrated using estimated en- 
ergy losses of e's, /i's and 7r's that passed through the 
counters. The telescope counters were calibrated with 
positrons from fi + — > e + uv decays, 7r+ — > e + u decays and 
scattered beam particles. For TINA and MINA, energy 
losses in the telescope counters and the 1-MeV shift due 
to positron annihilation were taken into account. Time 
calibration was based on the repeated peaks with the RF 
cycle in the time spectra. Figure 4 shows a scatter plot of 
energy in the T2 counter vs energy in TINA, in which the 
diagonal band is from protons produced by pion capture, 
the flat distribution around the energy loss of 2 MeV/cm 
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is from electrons (up to 45 MeV from /i — > evv decays and 
around 60 MeV from scattered beam electrons), and the 
distribution around 100 MeV with twice the energy loss 
of electrons in the T2 counter is from electron-positron 
pairs from 7-rays converted in the Tl counter or the tar- 
get. 

A. Analysis of -n — > /i — > e decay 

In the analysis of 7r~ — > fi~ — > e~ decay events, data 
taken with the circular water target were used. Pulse 
heights in the three telescope counters T1-T3 were re- 
quired to be less than ~7 MeV/cm. The effect of the 
loss by this cut (^3 %) cancelled when the yield was 
normalized to that of tt + — > fi + — > e + . No proton con- 
tamination below E=80 MeV was observed. Some events 
from two "electrons" were accepted with the dE/dx cut 
in the T-counters, but they appeared together with scat- 
tered beam electrons at the prompt time and at every 43 
ns with respect to the prompt time. Only delayed events 
outside those affected regions (i e = 48-81 ns with respect 
to the prompt time to)[l3 were selected for the analysis. 
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Energy in B2 (MeV) 

FIG. 5: Energy spectrum measured in the B2 counter. The 
dotted line shows the result of MC calculation normalized to 
the total count in the spectrum. 

The histogram of Figure 5 shows the energy spectrum 
measured in the B2 counter for the events that passed the 
above cuts. The position of the B2 counter was in contact 
with the target (z=0 cm). DIF contributions before and 
after passing the B2 counter are the major remaining 
components. Events from decays after passing the B2 
counter form a peak around 5 MeV that corresponds to 
energy loss of pions passing through the B2 counter prior 
to stopping in the target. Events from the decays before 



passing the B2 counter form a broad distribution with a 
shoulder around 3.5 MeV that corresponds to the energy 
loss of the highest energy muons (^32 MeV for forward 
pion decays) and the lower energy events are from the 
backward-decay muons that stopped in the B2 counter. 

The contribution of DIF events was simulated by MC 
based on GEANT3 13] with a modification for the 7r~ 
treatment. Pions were generated at the exit of the beam 
pipe (110 cm upstream of the target) with a profile of 
10 x 10 cm 2 projecting to the target (2x2 cm 2 ) accord- 
ing to a flat momentum distribution centered at 70 (± 
4 %) MeV/c that reproduced the measured energy spec- 
trum of pions taken with a beam trigger. When a tt~ 
stopped in material, it was promptly captured (capture 
products were not generated in the n — > // — > e simula- 
tion). Muon capture was not included in the MC, but the 
effect was included in the lifetime correction by using the 
free muon lifetime instead of the lifetime of the muonic 
atom The effect of the binding energy (0.2 MeV for 
the muonic Oxygen atom) on the electron spectrum was 
ignored. The MC events were analyzed with the same se- 
lection criteria as the experimental data and the energy 
spectrum in the B2 counter was obtained as shown by the 
dotted line in Fig. 5 (normalized to the total count). The 
shape of the distribution is well reproduced by MC cal- 
culations of the DIF events described above. The ratio of 
the two components in the MC spectrum, the peak (DIF 
after B2) and the broad bump (DIF before B2), strongly 
depended on the gap z, especially near the target (z—0 
cm). Although the peak is, in principle, more sensitive 
to the free decay component, the sum was used in the 
analysis because of this strong dependence near z—0 cm 
for reliable subtraction of the background estimated by 
the MC calculations. The contamination from the /z's in 
the beam was negligible. 

In order to obtain the free decay component fraction, 
the yield was normalized to that of the n + — > /i + — > e + 
decays in the same time window with the corrections for 
the lifetime and the stopping fractions obtained from the 
beam-trigger data (f v + : f v - = 0.81 : 0.69) 45]. The 
second column in Table I shows the gap dependence of 
the free decay fractions, which are dominated by the DIF 
contribution. For the normalization of the MC calcula- 
tions, 7r + data were also generated and analyzed in the 
same way. The last column in Table I summarizes the 
corresponding fractions from the MC calculations. 

The two measurements at z = cm represent data 
with a possible beam position shift due to the magnet 
setup adjustment to correct for the hysteresis effect. 
The average of two measurements at z=0 cm was used 
as the final result. An upper limit based on Bayesian 
method assuming a Gaussian probability distribution 
was obtained by subtracting the MC result and adding 
the associated errors in quadrature, and then renormal- 
izing the area above zero. The obtained 90 % Confidence 
Level (C.L.) upper limit of the n — > /1 — > e amplitude 
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Measured (%) 


MC (%) 




Material 


7r — ► ev 


7T Stops 


tv — > ev 


7T Stops 


90 % C.L. 





1.12 ± 0.04 


1.09 ± 0.03 




H 2 


3 


1.29 x 10 9 


69 


5.23 x 10 7 


4.6 x 10" 3 




1.09 ± 0.04 






Be 


1 


8.32 x 10 8 


94 


5.40 x 10 7 


3.1 x 10~ 3 


1.4 


1.26 ± 0.05 


1.17 ± 0.04 




Al 


1 


3.34 x 10 8 


114 


6.38 x 10 7 


7.6 x 10~ 3 


3.5 


1.15 ± 0.04 


1.07 ± 0.03 




Vinegar 





4.51 x 10 8 


69 


5.23 x 10 7 


4.5 x 10" 3 



TABLE I: Summary of n — > /i — ► e analysis for the water tar- 
get. To obtain the free decay fractions for the experimental 
data and the MC calculations, the yields have been normal- 
ized to those from tv + , and the time- window correction has 
been applied. 

corresponds to free-decay fraction of 7.5 x 1CT 4 for 
a short lifetime and 8.2 x 10~ 4 for a lifetime of the 
free-decay component near T b s = T n . 



TABLE II: Summary of tt -> ev analysis. The 90 % C.L. 
upper limits of the free decay fractions are shown in the last 
column for the case when t (, s = tv . 

II. The 90 % C.L. limits shown in the Table correspond 
to the case when the free decay component has the pion 
lifetime. These limits are complementary to the capture 
measurements. 



B. Analysis of tt — ► ev decay 

The signal region was defined by a box cut of i e =5-100 
ns and E e =60~72 MeV in TINA. Event selection criteria 
for the 7T — ► ev decay analysis were tighter than those 
used in the tt — > fj, — > e analysis because of its small 
branching ratio and sensitivity to background. A typical 
type of background occurred when a pion stopped in the 
target and an additional beam electron was scattered into 
TINA. Protection against this came from efficient beam 
pile-up detection and the high TINA energy threshold — 
the presence of extra beam particles during the search 
time window was detected by tighter energy loss cuts 
in the Bl and B2 counters as well as by comparing the 
pulse charges observed in the narrow-gate and wide-gate 
ADC's for the Bl counter. Using the observed pile-up 
data and beam electron energy spectra, the level of this 
background in terms of free decay fraction was estimated 
to be 4 x 1CT 3 . Pile-up in TINA that shifts the energy 
of tt — > fi — > e events into the signal region was another 
concern. Early pile-up in TINA was detected by the ADC 
with an early gate that closed at the beginning of the 
normal signal. Other pile-up was detected by a set of 
two ADC's with different gate widths for TINA and T3. 
This background was < 1 % of tt + — > e + u decays and 
significantly suppressed for n~ by the lack of free pion 
decays. 

Figure 6 shows energy spectra in TINA for tt + (top) 
and 7r~ (bottom) runs with the water target. Some 
events were seen in the time-energy box with time dis- 
tribution consistent with the electron scattering back- 
ground, though they were included in the free-decay can- 
didates. 90 % C.L. upper limits were obtained assuming 
Poisson statistics with no background. The limits were 
normalized to the corresponding tt + data with the time 
window corrections. The results are summarized in Table 




TINA Energy (MeV) 

FIG. 6: Energy spectra in TINA for tt + (top) and 7r~ (bot- 
tom) stopping in H2O. Events in the energy region of 60-70 
MeV are mainly from tt — > ev decays. 



C. Analysis of capture products 

In the search for free tt~ decay via detection of delayed 
protons from tt" capture, good timing is essential to be 
sensitive to fast decay components. In order to mini- 
mize the time smearing, the energy loss as well as the 
TOF cuts in the selection of pions (typically ±3cr away 
from the peak) were tighter than those used in the elec- 
tron analysis. Selection of protons was made based on 
the energy loss in all the three T counters; the thresh- 
old for each counter was set around four times of that of 
minimum ionizing particles. The resulting cut-off energy 
for accepted protons was around 45 MeV in the incident 
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energy. The contributions from electron and muon con- 
taminations in the beam were negligible. 
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FIG. 7: Time spectra of protons for H2O, Be, Al and vinegar 
(histograms) . The solid curves show Gaussian fits. 

Without corrections the time resolution of the prompt 
peak was a ~1 ns. There were many effects that 
caused timing shift and spread; 1) the beam momentum 
distribution that affected the energy loss in the beam 
counters, the horizontal stopping position, and the TOF, 
2) the pulse height dependences (time walk) in the beam 
counters and telescope counters, which were correlated 
with the energy loss and the proton emission angle, 3) 
the path length variation, e.g. due to the pion stopping 
position and the emission angle of protons, and 4) the 
propagation time of the light in the scintillator especially 
in telescope counters, Tl, T2, and T3. Many effects were 
correlated with one another. Correction factors were 
obtained by minimizing the dependence of the in on each 
measured parameter, such as pulse heights, timing in the 
beam counters and telescope counters, and horizontal 
pion stopping positions. A linear dependence was used 
to obtain corrections. New correction factors were 
calculated at the end of each iteration and the revised 
correction factors were applied in the next iteration 
(only half of the amount was corrected each time). 
The process converged after about ten iterations. The 
corrections improved the timing resolution to a ~ 200 
ps. These corrections do not affect the amount of the 
real decay component since the free decay component 
fraction is independent of the parameters used in the 
optimization. 

Figure 7 shows time spectra for H2O, Beryllium, 




FIG. 8: -Log(Lk) for H2O and Beryllium. Large — logLfc 
values were truncated for display. 



Aluminum and vinegar. The time spectra of protons 
for all targets show essentially a single component, a 
Gaussian peak at in 7 & s indicated by the solid line for a 
single Gaussian fit with three parameters; the amplitude, 
offset and width, reflecting the resolution. In the data 
from early stage of the experiment, there was an obvious 
delayed component, a bump at t p =l 2 ns that was 
enhanced when the events stopping in the upper half of 
the target were selected. The appearance of the delayed 
component coincided with the runs with the water target 
filled only up to 85 % of the height. These suggested 
some slow pions stopped outside the target (probably 
downstream of the target) with a longer flight path for 
the resulting proton. MC calculations indicated that, 
when the target was lowered by the half height of the 
beam, about 0.3 % pions stopped on the surface of the 
material around the telescope counters and had a similar 
time distribution. Since the extra counts outside the 
Gaussian distribution could not positively be identified 
as background, we treated them as a potential signal 
and only upper limits are quoted. 
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FIG. 9: 90 % C.L. upper limits of delayed capture ratios f cap t 
for H2O and Beryllium as a function of assumed lifetime. The 
true 90 % C.L. limits around the 1-ns region may lie along a 
curve that smoothly connects the binning and fit results. 

For the extraction of a possible delayed component by 
the time spectrum fit, the function used in the fit was a 
Gaussian plus an exponential function including the res- 
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olution effect (the same resolution was used). Assuming 
only one decay constant for the delayed component, a 
likelihood function Lk 



Lk 



n 



N n t 



-Ni n 

-/n 



(3) 



Material 


Candidates 


Total events 


90 % C.L. ratios 


H 2 





1.31 x 10 5 


2.3 x 10" 5 


Be 


1 


7.35 x 10 4 


6.9 x 10~ 5 


Al 





4.08 x 10 4 


7.3 x 10~ 5 


Vinegar 


1 


6.77 x 10 4 


7.5 x 10" 5 



where n% and Ni are the data and fitted value for the 
i-th bin, respectively, was maximized with three variable 
parameters; the total counts, the resolution, and the to, 
but the lifetime and the capture ratio of the delayed 
component were fixed for each maximization process, 
and the minimum in — log Lk was stored for each set of 
lifetime and delayed capture ratio. Fig. 8 shows surface 
plots of the likelihood functions — log Lk for the delayed 
capture ratio and the lifetime. The dotted histograms in 
Fig. 9 show the 90 % C.L. upper limits of the delayed 
capture ratio f ca pt for each lifetime. 

In order to test the validity of using a single Gaussian 
distribution for the prompt capture timing peak, a possi- 
ble asymmetric time distribution, due to e.g. solid angle 
effects causing different weight for different TOF, was 
simulated by MC. Capture products were generated as 
soon as they stopped; pion capture in flight was simulated 
separately. Proton and deuteron spectra were taken from 
Ref.|l6l] and interpolated in energy and atomic number. 
Other heavier particles were ignored because of higher 
stopping power that prevented them to contribute to a 
trigger. The MC events were calibrated in the same way 
as the experimental data to reflect the effects of the cor- 
rections in timing spectrum, and analyzed with the same 
program with the same cuts. After smearing with a reso- 
lution of 200 ps (the bare resolution of MC was about 60 
ps after applying corrections), the spectrum was fitted to 
a Gaussian curve. The fit agreed at the 10 _4 -10 -5 level, 
and the asymmetry was not significant. 

The effect of capture-in-flight was simulated with 
the total cross section 0, a tota i = Trrg ' (1 + "f^")' 
where ro = 1.25 • A 1 / 3 (/m) is the nuclear radius, and 
E c = 1-44 • Z/r (MeV) is the Coulomb potential. A 
and Z are the mass and atomic numbers of the target 
material (Hydrogen was ignored for scintillator and wa- 
ter), respectively. The nuclear absorption cross section 
was assumed to be 1/3 of the total cross section, and the 
capture products were generated by ignoring the kinetic 
energy of the pion. Inclusion of the capture-in-flight 
process increased the number of events in the early time 
region (t p = —0.5 ns) by 0.1 % of the total counts, which 
was consistent with the excess counts in the same time 
region of the data. However, the capture ratio of the 
delayed component was unaffected by this. 

In the long lifetime region, the fit results were affected 
by the small sample bias of the Maximum Likelihood 
Method. A simple time bin method was employed for 



TABLE III: Summary of binning analysis of proton data. The 
second column shows the numbers of candidate events in the 
time bin 5-75 ns. The 90 % C.L. upper limits for delayed 
capture ratios in the last column are for r bs = tv with the 
corresponding time- window correction of 0.769. 



this region; events in the time bin of t p = 5 — 75 ns 
were considered to be possible candidates, and the sum 
in the time bin was normalized to the total number of 
events in the proton spectrum. A Poisson distribution 
was assumed to obtain the 90 % C.L. upper limit. Table 
III summarizes the results for H2O, Be, Al and vinegar. 
The solid lines in Fig. 9 show the results from the binning 
method. 
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FIG. 10: 90 % C.L. upper limits of free decay fractions for 
H2O, Be, Al, and vinegar. The dashed histograms are the 
results by fit, the solid histograms by binning for the delayed 
capture method, the dotted histograms by n — > ev and the 
dash-dotted line by 7r — > ji — > e. 

The observed delayed-capture ratio was converted to 
the fraction of free ir~ decay using Eq.(2). The summary 
of the proton analysis is shown in Fig. 10, where the 
dashed lines indicate the 90 % C.L. limits for the free 
decay component obtained by the fits, and the solid lines 
by the binning method. The capture limits go to unity 



8 



at the pion lifetime (the next bin outside the spectrum). 
In the shorter lifetime region outside the plot, the upper 
limit is along the extension of the binning result, which 
corresponds to r hal T -n- The results from the n — ► )jl — ► e 
(dash-dotted line) and n — ► ev (dotted line) analyses are 
also shown in the figure. Without knowing the lifetime 
of the delayed component, conservative limits obtained 
near the pion lifetime (the largest upper limits) for all 
cases are: 8.2 x 10~ 4 for H2O by the tt — > fi — > e analysis 
at r obs = t„, and 3.2 x 10~ 3 for Beryllium, 7.7 x 1(T 3 
for Aluminum and 4.5 x 10~ 3 for vinegar by tt — > ev at 
the intersection of the binning and tt — > ev limits. 



V. CONCLUSION 

No free 7r~-decay component was observed in water 
with the "90 % C.L." upper limit for the free decay 
fraction of 8.2 x 10~ 4 for the worst case (t i, s = r n ). 
For Beryllium and Aluminum assuming a single delayed 
component, the limits were 3.2 x 10~ 3 and 7.7 x 10~ 3 
coming from the measurements of 7r~ — > e~v decays, 
respectively. The result of the present experiment indi- 
cates that the materials studied here will not cause v e 
contamination in the beam dump neutrino experiments. 
It seems unlikely that the LSND result 0, which would 
require the relative v e yield of > 10~ 2 ^3j can be 
explained by the mechanism of unusual formation or 
de-excitation processes of pionic atoms in water. 
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